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ABSTRACT 

The objective of this paper is to present a comprehensive approach to 
analyzing the dynamic behavior of multistage gear transmission systems, with 
the effects of gear-box-induced vibrations and rotor mass-imbalances. The 
model method, using undamped frequencies and planar mode shapes, is used to 
reduce the degree-of-f reedom of the system. The various rotor-bearing stages 
as well as the lateral and torsional vibrations of each individual stage are 
coupled through localized gear-mesh-tooth interactions. Gear-box vibrations 
are coupled to the gear stage dynamics through bearing support forces. 
Transient and steady state dynamics of lateral and torsional vibrations of the 
geared system are examined in both time and frequency domains. Vibration 
signature analysis techniques will be developed in interpret the overall 
system dynamics and individual modal excitations under various operating 
conditions. A typical three-staged geared system is used as an example. 
Effects of mass-imbalance and gear box vibrations on the system dynamic 
behavior are presented in terms of modal excitation functions for both lateral 
and torsional vibrations. Operational characteristics and conclusions are 


drawn from the results presented. 
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INTRODUCTION 


Presently there is a continuous search for the improvement in 
operational life, efficiency, maintainability, and reliability in gear 
transmission systems. With the increase in workload and speed requirements, 
the call for transmission design improvement becomes even greater. One of the 
major objectives in design improvement is the reduction of noise and vibration 
in the transmission system. Two main streams of work have been carried out in 
the areas of noise and vibration reduction; namely, (1) the localized gear 
tooth stress and thermal effects during gear interactions (Cornell 1981; Lin 
1988; Boyd 1987; Savage 1986) and (2) the overall global dynamic behavior 
(August 1986; Choy 1989? David 1987, 1988; Mitchell 1985) of the transmission 
systems . 

The work presented in this paper is the development and application of a 
comprehensive approach in simulating the overall dynamics of a multistage 
rotor-bearing-gear system. A system of rotor-bearing configurations with 
circular flexible shafts on flexible bearing supports are considered in this 
analysis. Excitations due to rotor mass-imbalance, shaft bow, and gyroscopic 
effects are considered in each individual rotor stage. Effects of base motion 
are incorporated through relative motion between rotor and casing at the 
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bearing supports. The dynamics of each rotor stage are coupled together at 
gear locations through nonlinear spring connections. The modal method is 
applied to reduce the degrees-of-f reedom of the global system of equations 
into modal coordinates. The modal equations of motion are solved to evaluate 
system acceleration at each time step. A self-adaptive variable time-stepping 
integration technique (Choy 1987, 1988a, 1989) is used to calculate the 
transient dynamics of the system. A typical three-stage rotor-bearing-gear 
system is used as an example in this analysis. Results are presented in both 
time and frequency domains to develop vibration signature analysis of the 
system. 

DEVELOPMENT OF EQUATIONS OF MOTION 

The equations of motion for a single stage multimass rotor-bearing-gear 
system with the effects of gear-box-motion-induced vibrations at the bearing 
supports, rotor inert ia-gyroscopic effects, and excitations from rotor mass- 
imbalance and residual runouts can be written in matrix form (Choy 1987, 1989) 
for the X-Z plane as 

[M](x) + [G v ](i} + [C xx ](* - * b } + (C xy ]{i - t b ) + (G a ]{y) 

+ l K xx + M*) - f K xxj( X b) + l K xy) ( Y - Y b} = {*.<*>} + { F GX ( fc > } 

and in the Y-Z plane as 

[M](?} - [o v j{*} + [c yx ]|X - * b ) + [c yy ](* - i b } - [g a ](x} 

( 2 ) 

+ [Kyy + K J{ Y ) - t K yyl { Y b} + t R yx] X ( Y " Y b) " { F y < > } + { F Gy < > } 

Here F and F are force excitations from the effects of mass-imbalance 

x y 

and shaft residual bow in both x and y directions. Forces F qx and are 
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the X and Y gear mesh forces induced from the gear teeth interaction with 
other coupled gear stages. The bearing forces are evaluated through the 
relative motion between the rotor {X}, {Y} and the gear box {X^} , { Y^} at the 

bearing locations (Choy 1987), The mass-inertia and gyroscopic effects are 
incorporated in the mass matrix [M] and the gyroscopic matrices [G v ] and [ G a } - 
The coupled torsional equations of motion for the single rotor-bearing-gear 
system can be written as 

[J]{e) + [C T ]{^} + [K t ]{0} = {f t ( t ) } +■ (F Gt (t)} < 3 > 

In Eq. (3) {F (t)} represents the externally applied torque, and {F Qt (t)} 
represents the gear mesh induced moment. Note that Eqs. (1) to (3) repeat for 
each stage. The gear mesh forces couple the force equations of each stage to 
each other as well as the torsional equations to the lateral equations 
(Choy 1989; Cornell 1981; David 1987, 1988). The coupling relationships 
between the torsional and the lateral vibrations and the dynamics of each 
individual gear/rotor are derived in the next section. 

COUPLING IN GEAR MESHES 

The torsional and lateral vibration of each individual rotor stage and 
the dynamic relationships between all the gear stages are coupled through the 
nonlinear interactions in the gear mesh. Gear mesh forces and moments are 
evaluated as functions of relative motion and rotation between two meshing 
gears and the corresponding gear mesh stiffnesses. These gear mesh 
stiffnesses vary in a repeating nonlinear pattern (August 1986; Cornell 1981; 
Savage 1986) with each tooth pass engagement period and can be represented by 
a high order polynomial (Cornell 1981; Boyd 1987). A sixth order polynomial 
compliance curve is used in this study to simulate the stiffness changes for 
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contacting gear pairs (zero stiffness is input for noncontacting pairs). The 
repeat ibility of such nonlinear mesh stiffnesses can also act as a source of 
steady state excitation to the gear system. With the coordinate system as 
shown in Fig. 1 and with the effects of gear tooth surface friction neglected, 
the following gear mesh coupling equations can be established by equating 
force and moment in terms of base circle radius R (Choy 1989). For the k th 

c 

stage gear of the system, summing force in the x direction results in 


? G xk = E K tlJ- R cAi - R cAk + < X ci - X ck> COS «ki + < Y cl - Y c k ) Sin a Ki] COS Q ki 

i=I,ii*k (4) 


Summing force in the y direction results in 
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i=l r ij*k (t 
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( 6 ) 

where n is the number of stages in the system. Damping forces in a similar 
form involving the relative velocities between the two mating gear stages 
would also be included. 

MODAL ANALYSIS 

To reduce the computational effort, the number of degrees-of-f reedom of 
the system is reduced through modal transformation. Orthonormal modes for 
each individual rotor-bearing stage are obtained by solving the uncoupled 
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system homogeneous characteristic equations. Using the modal expansion 
approach (Choy 1987, 1988a, 1989), the motion of the system can be expressed 


as 


M ■ E 
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where m is the number of modes used to define each motion, 
orthogonality conditions of the modes can be expressed as 


[$] t [K][$] = tA 2 ] 


where 
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( 8 ) 
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[* t J T [K T ][* t ] = [Aj 


( 11 ) 


[$] t [M][$] = [$ t ] T [J)[* t J = [I] 


( 12 ) 


Using the modal expansion and the orthogonality conditions, with the bearing 
forces due the base motion expressed in the R.H.S., the modal equations of 
motion (Choy 1989) can be written: 

For the X-Z equation 


(a) + [$) T [G v ][$]{s} ♦ (*] t [c xx ]($](a) + [* ] T [c xy ][*](6} + [*) t [g a ]($]{b} + [A 2 j ( a) 

+ [*] T [K dx ][$](A} + [*] T (K xy ][*]{B} = {*] T {F x (t> +F 0 x (t) +F ex (t)} 

(13) 


where 

F B X < fc ) = [ C xj(*b) + [ C xy)( Y b} + f K xx>( X b} + f K xyj( Y b}( K dxJ = ( K xx 1 “ l K l 


For the Y-Z equation 
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where 


F sy (t) = [^{Xj - lC yy ]{Y t ] + [ K yx ](X b ) * [K yy \{Y b )[K dy \ = [K^] - [K] 


And for the 0-equation 
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( 15 ) 


The gear mesh forces and moments can also be expressed in the modal form, for 
the k th stage with gear location at the 1 th node, as 
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where k is the stage number, j is the mode number, and 1 is the station 
location of the gear mesh. 

SOLUTION PROCEDURE 

Rearrange the modal equations of motion developed in Eqs. (13) to (15) 
into the following forms: 

X-equat ion 


(a) = -1$] t [G v h$](b} - [$] t [C xx ][$j(a) - [ $ ] T [ c xy ] [ $ ] { B } - [$J T [G a ][$](b} - [A 2 ] {a} 
- [$] T [K dx ] [*]{a) - [$] T [K x ..][$]{b} + [$] T (F x (t) +F Gx (t) +F Bx {t)j 

(19) 


Y-equat ion 


(b) = [$] t [g v ][$)(a) - ($] t [c yx ][$)(a) - [*] T [c yy ][*]W + [$] t [g a ][$]{a) - [A 2 ){bJ 

X t^] T + [ K dy J I $ ] ( B ) - [•] T [K yx H* ] {A} + [$] 7 (F y (t) + F Gy (t) + F By (t)} 

( 20 ) 


0-equation 


(A t ) = -[$ t ] T [C T ] [* t ](A t } - [Aj(A t } + ($ t ]’ I (F t (t) + F Gt (t)j 


( 21 ) 


A variable time stepping Newmark— Beta integration scheme evaluates the modal 
velocity and displacement at each time interval for each stage. In turn, 

Eqs. (7) to (9) transform the modal displacements into absolute displacements 
in fixed coordinates. The gear mesh forces can be evaluated by the relative 
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motion between the gear teeth using the nonlinear stiffnesses developed for 
the gear mesh interaction. The effects of gear box vibration can be 
calculated as nonlinear bearing forces through the relative vibration between 
the gear box and the shaft at the bearing locations. Since the gear box and 
the shaft are vibrating independently, a separate transient integration scheme 
is required for each system at every time step before the coupling of bearing 
forces . 

DISCUSSION OF RESULTS 

To demonstrate the application of the discussed analytical method, a 
typical three-stage gear transmission, given in Fig. 2, is used as an example. 
All three gear stages have an identical 36 tooth gear and a mesh contact ratio 
of 1.6. Stage I of the system is the driver with a rotational speed of 
1500 rpm and an input torque of 2.25 kN*m. Both stages I and II are supported 
by three bearings; stage III only has two bearing supports. The first two 
bearings in stages I to III and the third bearing in stages I and II are 
identical. The rotors in stages I and II are of 2.4 in. in diameter and are 
of length of 40 and 30 in., respectively. The rotor in stage III is smaller 
with a diameter of 1.4 in. and a length of 20 in. For demonstration purposes, 
the first three basic lateral and two torsional modes of the system are used 
in the modal anlaysis (numerical experiment shows that the first three basic 
lateral modes accounts for over 90 percent of the system vibrations) . Table I 
provides the information on the modal frequencies of the system in both 
lateral and torsional directions. Figures 3 and 4 give the orthonormal mode 
shapes of the first two basic torsional and lateral modes for all three rotor 
stages. (Third lateral mode shape and its modal excitations are omitted from 
the figures because of the mode shape's small magnitude in vibration). 


11 



In order to investigate the effects of gear box vibrations and rotor 
mass-imbalance, results from four major cases of external excitations are 
examined in this study: 

Case 1: no gear box vibration and small mass-imbalance 

Case 2: gear box vibration in x-direction and small mass-imbalance 

Case 3: no gear box vibration and large mass-imbalance 

Case 4: gear box vibration in x-direction and large mass-imbalance 

The gear box vibrational effects are assumed to be a steady state 
vibratory motion of 600 Hz in the x direction at the bearing supports and to 
be unaffected by shaft vibrations. In addition to the damping characteristics 
due to the bearing supports, a 1-percent modal damping is assumed at the gear 
mesh. An initial condition of 0 velocity and displacement is used in this 
study. Vibrational orbits of all three rotor stages at the gear location for 
all four cases are given in Figs. 5 to 7 . Note that the vibratory 
characteristics of stage I are more influenced by mass-imbalance while those 
of stages II and III are more affected by motion of the gear box thorugh the 
bearing supports. As we can see from case 2 of Fig. 5{b), the effects of 
bearing support motion on stage I gear vibrational orbits are quite small. 

This is due to the fact that the rotor of stage I is substantially more 
flexible than the other two rotors and less influence from the bearing can be 
effectively converted to the system. In addition, although similar mass- 
imbalances are applied to each rotor stage, vibratory orbits in stages II and 
III are smaller in amplitude because of their larger rotor stiffnesses. 

Figures 8 to 10 depict the x direction first modal component for rotor 
stages I to III, respectively. Figure 8 shows the effects of the four 
external excitations on the first modal vibratory component of stage I in both 
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time and frequency domains. Note that Fig. 8(a) shows that, for small 
imbalance with no gear box motion, only very small vibratory motion will 
result, except for the dc component at 0 frequency (due to the x direction 
displacement of the center of the orbit from the original journal position). 
Figures 8(b) represent the first mode rotor vibrations resulting from gear box 
motion. Note that the major component occurs at 600 Hz (input vibratory 
frequency at the bearing supports) with a minor component at 115 Hz (first 
natural frequency of stage I). Figures 8(c) show the effect of mass-imbalance 
which results in a large component at the rotor speed of 25 Hz and a smaller 
component at the first mode of 115 Hz. Figures 8(d) produce the combined 
effect of gear-box vibration and imbalance on rotor first mode vibrations. By 
examining the relative amplitudes of the two major components, at rotor speed 
and input frequency, with their corresponding input excitation magnitude, the 
sensitivity of the rotor to imbalance and gear-box vibration can be 
determined. Similar conclusions can be arrived in Figs. 9 and 10 for stages 
II and III except that the effects of mass-imbalance on first mode vibration 
is significantly reduced. As we can see in Figs. 9(d) and 10(d), when both 
excitations are presented, the gear-box motion has a much more dominating 
effect in rotor vibrations for stages II and III then for stage I. 

Figures 11 to 13 show the vibration of the second mode component for the 
three rotor stages. Note that for stage I, the second mode vibrational 
characteristics are very similar to those of the first mode (Fig. 11). 

However, actual magnitudes are approximately 20 percent lower than those of 
the first mode. This is due to the fact that the first mode is more easily 
excited by the imbalance as well as by the bearing support excitations. 

Similar trends of vibrational effects can also be seen in the rotor motions in 
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stages II and III (Figs, 12 and 13), Note that the second modal frequency is 
not excited in this case, which further confirms the high excitability of the 
first mode. In addition, none of the torsional vibratory frequencies and gear 
meshing frequencies are excited in the lateral rotor vibration of the system 
( Figs . Q to 13 ) . 

Figure 14 depicts the gear mesh forces with effects of both imbalance 
and gear-box motion excitations. Note that the major component of these 
forces are the dc component at 0 frequency due to the constant applying 
torque. The other sizable force component is at the tooth meshing frequency 
of 900 Hz. The higher ac component from the meshing of gears I and III 
results from the lower rotational stiffness in the stage III rotor. The 
effect of this lower rotational stiffness is also shown by the higher dc 
(0 frequency) component in stage III compared with stage II (Figs. 15 and 16). 
Only very small trace of force components can be observed at the gear mesh 
frequency (900 Hz) and first mode torsional frequencies of stage I (355 Hz), 
stage II (550 Hz), and stage III (280 Hz). Again, Figs. 15 and 16 outline the 
modal rotational vibration characteristics of all three rotor stages with both 
excitations. Note in Fig. 15 that other than the 0 frequency component, the 
other sizable vibration is at its own rotational modal frequency (i.e., 355 Hz 
at stage I, 550 Hz at stage II, and 280 Hz at stage III). Figure 16 depicts 
the second mode rotational vibration characteristics. Note that similar 
conclusions can be reached for the second modal frequencies of 1090, 1610, and 
820 Hz, respectively, for stages I to III. 

SUMMARY 

This paper presents a vibration analysis with the effects of gear box 
motion and mass-imbalance for a multistage gear transmission. The analysis 
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combines gear mesh dynamics and structural modal analysis to study the 
transmission vibrations. The major content of this work can be summarized as 
follows : 

1. A comprehensive approach is developed to combine the nonlinear gear 
mesh dynamics with structural lateral and torsional vibration of the system to 
determine the global system response. 

2. The modal method transforms the equations of motion into modal 
coordinates to reduce the degree-of-f reedom of the system. 

3. Gear force observations in both the time and frequency domains 
provide good insights into the source of dominating response forces. 

4. The magnitude of the gear mesh force is inversely proportional to the 
rotor stiffness of the driven system. 

5. The influence of the gear box motion on system vibration is more 
pronounced in a stiffer rotor system. 

6. Gear tooth mesh frequency and torsional modal frequencies have 
substantial effects on rotational but not on lateral vibrations of the system. 

7. Knowledge of modal amplifications under various excitations provide 
an understanding of the vibrational character istics of the system. This 
knowledge is crucial for designing transmissions with improved performance and 
durability. 
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TABLE I. - SYSTEM NATURAL 
FREQUENCIES 


Mode 

Stage 1 

Stage 2 

Stage 3 

Torsion natural frequency 

1 



280 

2 



820 

Lateral natural frequency, Hz 

1 

mm 

160 

110 

2 


189 

200 

3 

190 

264 

260 
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Figure 4.— Second vibrational mode for the three-stage system. 
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Figure 5. — Stage 1 rotor orbits 
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Figure 7.— -Stage 3 rotor orbits under various load conditions. 
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Figure 14. — Gear forces in time and frequency domains. 
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dynamics through bearing support forces. Transient and steady state dynamics of lateral and torsional vibrations 
of the geared system are examined in both time and frequency domains. Vibration signature analysis techniques 
will be developed in interpret the overall systems dynamics and individual modal excitations under various 
operating conditions. A typical three-staged geared system is used as an example. Effects of mass- imbalance and 
gear box vibrations on the system dynamic behavior are presented in terms of modal excitation functions for both 
lateral and torsional vibrations. Operational characteristics and conclusions are drawn form the results presented. 
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